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Abstract 

We provide an alternative derivation of a lower bound on the mass of the Higgs boson 
which is somewhat simpler and more direct than the derivation based on the effective 
potential. For one TeV cutoff, the result is the same. For high scale cutoff, the lower 
bound is increased by slightly more than the expected uncertainty in the calculation. 
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For the history of the subject we refer to reviews [[TJ and quote only two recent papers 
containing the latest refinements of the conventional approach ||. Once the experimental 



lower bound on the top-quark mass exceeded 80 GeV, attention shifted from the original 
Linde- Weinberg bound, based on the properties of the one-loop effective potential for 
small <p, close to the minimum, to the large behavior of the effective potential , as 
determined by renormalization group (RG) considerations. 



Here, is the anomalous dimension factor, and A(t) is the MS running coupling con- 
stant . 



It is then argued that vacuum stability requires A (t) > up to some high scale, 4> ~ 
Mqut, or Mpi, (M ~ Mz, or m t , or 246GeV).To implement this condition, one has to 
know (or approximate) the (3— function, integrate the RG differential equations starting 
from some initial values, A (0),g 2 (0), and relate the smallest acceptable A (0) to a physical 
Higgs mass. 

We do appreciate this calculation, but questions may be raised about the perturbative 
nature, the scale ambiguity, and the conceptual basis. The exact (3— function is not known, 
so one integrates the one- loop (3— functions to get the "RG improved one- loop effective 
potential " . Since the coupling is not asymptotically free, the large t behavior of A (t) is 
not known. The best that can be (and is) achieved is perturbative self-consistency. For 
the indicated A (0), integrating the one- loop (3— function gives A (t) which remains per- 
turbative up to the high scales considered. The minimum A (0) depends on the minimum 
scale Mo above. The physical Higgs mass does not. So there is some scale ambiguity. On 
the third point, we largely repeat the remarks of ||. If one thinks of a nonperturbative 
formulation of the vacuum stability problem, in particular, a lattice formulation; the large 
behavior of the effective potential is not the point. On the lattice, the exact effective 
potential is well defined and convex. The condition for vacuum stability is simply that 
the bare quartic coupling constant must be positive (Xbare > 0). Then the lower bound on 




dX(t) 
dt 
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the Higgs mass is just the smallest output Higgs mass from the Monte Carlo simulation 
as one runs through the space of bare parameters in the broken symmetry phase. 

We present a new derivation of the Higgs mass lower bound. It is also of the perturba- 
tive RG variety, and so also subject to the first concern above; but we have organized the 
calculation in a way which minimizes the scale ambiguity and makes no explicit reference 
to the effective potential . The essential input is that one is perturbing about the correct 
vacuum. A neccessary condition for this is that the vev of the (shifted) field be zero, order 
by order in perturbation theory, and the renormalized mass squared of the shifted field 
be positive. 

We start by computing the relation between the perturbative pole mass and the MS 
mass, for both the Higgs boson and the t-quark. The relation follows from the perurbative 
definition of the pole mass, 

= ~D~ l (M* 2 ) = M* 2 — ~M 2 — ReY,(M* 2 ) (1) 

In this equation, D(q 2 ), and S(g 2 ) are the two- point Green Function and self-energy 
function, renormalized according to the MS prescription. M* is the perturbative pole 
mass. The result is 



M* =M 2 {1+A(3/ 00 (M* )+9I jm (M* ))+N c f(^-A^)I Jfm (M^)+2(C v -l)} (2) 

m is the t-quark mass, and y is the t-quark Yukawa coupling (m = ^=). The contributions 
from the electroweak gauge sector, proportional to #2, have also been calculated, but 
are not written out here, they will be included below. The term ( v — 1 comes from a 
finite shift of the vev required in the MS scheme to enforce < H >= through one- loop 
order. It will cancel out of the ratio computed below, so we do not have to give its value 
here. ^I a b is the dimensionally regularized MS scalar one-loop two-point integral. 

T ( 2\ _ r A-d • r 1 1 
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Then (ffl) is 



M* 2 = M 2 {l + I ^[121n— -24+3V37r]+iV CI |^(l--)[ln-+/(r)]+2(C-l)} (4) 

3 



where 

M „, . U-r 2 I r 2 

r = — , /(r) = -2 + 2W — arctan \ / 

m V f y 4 — r z 

The corresponding calculation for the t-quark gives Q 

™* 2 = - 2 {! + | 1 n ^ + A W] + yI^HS - 6 In ^) + 2(C - 1)} (5) 

where 



A(r) = -4 + — + (-r 2 - ^r 4 ) lnr 2 + -(4 - r 2 )a arctan - ' 



2 ■ s 2 - 4' ' — ■ ' 2 V ~ ' V r 2 

We take the ratio of (0) to (H) and expand to one-loop order. 



M* 2 M 2 



■{1 + ^[12 In f - 24 + 3V3n] + ^[N c (l - £)(ln ^ + /(r)) - | In ^ - A(r)] 
+^C F (6 In s£ - 8) + «? 2 , ^ 2 terms + 2 - /oop} 

(6) 

The ( v — 1 terms, which also contain explicit dependence on ln/z 2 , have cancelled out. 
A necessary condition for the MS perturbation calculations to be defined in the broken 
symmetry phase is that M 2 ,rh 2 be positive. Since the ratio of pole masses is positive, 
(^D satisfies the requirement perturbatively, for fi around the weak scale. For large /tt 2 , 
one has to provide a RG treatment of the large logarithms, just as in the conventional 
calculation involving the effective potential . In the broken symmetry phase, one can 
define the renormalized coupling constants such that the relation 

K - 4 m 

m z y z 

is exact when all the quantities are either "star" (on-shell renormalization scheme) or 
"bar" (MS renormalization scheme) ||. Thus, not all quantities in (||) can be varied 
independently as functions of \x. We focus particularly on the multiplying the ln 2 ^. 
Tracing its origin to a ratio of MS masses in (||), we use (0) to replace the ratio of 
MS masses by a ratio of MS coupling constants. To leading (one-loop) order, the scale 
dependence of the ratio of MS masses is determined by the coefficients of the explicit 
In yU 2 terms in (Q). For the other masses in ([]), the difference between "star" and "bar" is 
higher order (combined with explicitly two-loop effects), as is the implicit \i dependence of 
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the "bar" coupling constants. After these observations, and reinstating the g%,gf terms, 
differentiating (|6]), we obtain 

f£® = it [ 24 (f) 2 + (2^c - 3 + 12CVff)|r - 2N C 

We now give a sequence of estimates, of increasing refinement, of the ratio Let 
=2 = p. (By (0), p = Let the right hand side of (§) be denoted /3 p . Because of 
the —2N C term in (0), there is a critical value of p below which (3 P becomes negative. 
And if the starting value of p is below this value, as p decreases the derivative becomes 
more negative, driving p negative, unless some higher order effect intervenes. We will 
return to this possibility, but as our zeroth order estimate we take the critical value of 
p for which one-loop f3 p , evaluated with weak scale coupling constants, is zero. We take 
g 2 s = 1.366 (a fl (174) = .109) and y 2 = 0.9990 [m t = 174, v = 246.2) and neglect the 
gig 2 contributions. Then p c = 0.2019, which gives (|J) 2 = 0.8076, or M c = 156, for 
m = 174. If we include the contribution from the electroweak gauge couplings, g%,gf, the 
corresponding results are p c = 0.2068, and M c = 158, a one percent shift. 

The first refinement is to convert back from (%r) c to the ratio of squared perturbative 
pole masses by (^J). Note that precisely for p = p c , all of the ln/x 2 terms in (^) cancel, so 
there is no explicit dependence on p in this correction. The result is 

M* 2 M 2 

Rr.= &c 1-0.101 (9) 

m* m 

which gives > 148. 

We now turn to the effect of the running of the MS coupling constants, which appear 
as coefficients in (Q),and the dependence of the lower bound on the cutoff (maximum 
value of -^-). One has to integrate coupled RG equations for five independent "coupling 
constants", g 2 s ,gl } g\ } y 2 ,p. Let t = ln£. 
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Itt- -3^(22 -\N a )t. 
m = ^[fN g + \N d ]g\ 

it = ^[(3 + 2 Nc )t - i2c F ttf - Ut - Mv 2 } 

|(p) = Xr [24p 2 + (2iV c - 3 + 12CV§|)p - 2iV c 

The first three equations are integrated trivially. If we neglect the g~2,gi contributions 
to the y running, that equation can also be integrated analytically. But if one runs 
up to high scales, the electroweak gauge couplings become of same order as the QCD 
coupling constant; so we use NDSolve from Mathematica to provide an interpolating 
function solution for y 2 which is subtituted into the g 2 equation, which is again integrated 
numerically by NDSolve. 

If we run up to the Planck scale (m p i ~ 10 19 GeV, t ~ 41) the smallest starting p(0) 
which does not lead to p(t) falling through zero before t = 41 is 0.2022, which is only 
slightly different from the critcal value required to make the derivative zero at the weak 
scale. At the other extreme, if we only require the equations of the standard model to be 
consistent up to order of one TeV, the value of t which corresponds to one Tev depends 
on the choice of [Aq. We choose liq = m t as the most natural choice for relating m to m*. 
For m = 174, this corresponds to t max = 1.75. Then the smallest starting p(0) which does 
not lead to p(t) falling through zero before t = 1.75 is 0.0520. Making the connection 
back to the ratio of pole masses by (Bl), we obtain the final result of this approach 

m* t = 174, a s {m t ) = .109 
AW « mpi M h > 148 (141, 135) (11) 

p max w ITeV M h > 72 (72) 

The numbers in parentheses are the corresponding results in the conventional aproach 
. The lower bound obtained in the present approach is slightly higher in the high scale 
cutoff case, but not by much more than the difference between the results of two different 
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calculations in the conventionl approach. 

The present derivation has the advantage that the zeroth order approximation, the 
value of — obtained for the vanishing of the one-loop (3 P at the weak scale, differs by 
less than ten percent from the final value for the large cutoff limit. If one makes the 
corresponding zeroth order determination of the minimum M in the conventional approach 
by setting f3\ to zero at the weak scale, the result differs from the final large cutoff result 
by more than thirty percent. The significant difference is that (3 P contains the large QCD 
correction to m, while (3\ does not. It is added in later as a correction when one integrates 
the coupled RG equations. 

It is clearly desirable to have a large scale lattice simulation study of the combined 
Higgs-heavy quark-QCD sector. (Contributions from light quarks and electroweak gauge 
bosons are small, particularly if one doesn't run up to some very high scale). We note 
that a quenched approximation simulation is not adequate for this problem. The term 
in (U) which triggers the possible instability is the —2N C , clearly a contribution from an 
internal fermion closed loop. 
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